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Abstract. Recent progress in numerical relativity now allows computation 
of the binary black hole merger, whereas post-Newtonian and pcrturbative 
techniques can be used to model the inspiral and ringdown phases. So far, most 
gravitational-wave searches have made use of various post-Ncwtonian-inspired 
templates to search for signals arising from the coalescence of compact binary 
objects. Ajith et al have produced hybrid waveforms for non-spinning binary 
black-hole systems which include the three stages of the coalescence process, and 
constructed from them phenomenological templates which capture the features of 
these waveforms in a parametrized form. As a first step towards extending the 
present inspiral searches to higher-mass binary black-hole systems, we have used 
these phenomenological waveforms in a search for numerically-simulated signals 
injected into synthetic LIGO data as part of the NINJA project. 



PACS numbers: 04.80.Nn, 04.30.Db, 04.25.Nx, 04.25.de 

1. Introduction 

An international network of ground-based gravitational-wave (GW) detectors 
(LIGO pQ, Virgo [2], GEO600 [3]) operating roughly in the 10 1 - 10 4 Hz band has 
recently finished taking data at or close to design sensitivity, and the world-wide 
LIGO Scientific Collaboration (LSC) and Virgo Collaboration are currently involved in 
analyzing it. Furthermore, space-based instruments, such as the Laser Interferometer 
Space Antenna (LISA) [HH], are also expected to probe a significantly lower- frequency 
stretch of the GW spectrum, namely the 10~ 4 — 10 _1 Hz band. One of the most 
promising candidate sources for detection is a compact binary system, composed of 
either neutron stars or black holes. GW emission from such systems, in particular 
from those consisting of binary black holes (BBH), is expected to be detected by 
these instruments in the near future, giving rise to the new field of gravitational-wave 
astronomy. 

The detection of gravitational signals arising from the coalescence of a BBH 
system, as well as the estimation of its physical parameters, relies heavily on the 
ability to accurately model the waveforms emitted during the inspiral, merger and 
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ringdown stages of the process. Until recently, the usual approach to the problem 
consisted entirely of making use of analytical approximations to the full theory of 
general relativity. Thus, a post-Newtonian (PN) expansion would be used to solve 
Einstein equations in the initial inspiral part of the coalescence (see [B] and references 
therein), while perturbation theory would be applied to the final ringdown after the 
two black holes have merged together. The intermediate phase corresponding to the 
strongest gravitational signal remained unmodcllcd until the breakthroughs occurred 
in 2005 pj [9] , when several numerical- relativity (NR) groups finally succeeded in 
using their numerical codes to solve for the merger and ringdown phases of the BBH 
coalescence, calculating the gravitational waves associated to the process. Steady 
progress quickly followed those pioneering results, leading to a thriving era in the 
field. 

Among the most suitable data-analysis (DA) methods currently used to look for 
signals from compact binary coalescences in the output data of GW interferometers 
is the matched-filter search technique against a given class of template waveforms. So 
far, most current LSC searches for compact binaries in the LIGO/ Virgo data have 
made use of different flavours of analytical PN templates to search for gravitational 
waves PHI [Hj ■ Further efforts have been undertaken in the recent past to incorporate 
the full inspiral-mcrgcr-ringdown (IMR) waveforms with information from NR results 
in the different DA packages into the LSC Algorithm Library [TJ]. These include 
the EOBNR waveforms based on the effective-one-body approach [13] and the 
phenomenological waveforms resulting from hybrid PN-NR match [TU [15j [16] among 
others. 

In order to be able to make confident statements about the performance of diverse 
DA pipelines and their ability to find true gravitational waves, data analysts have 
developed several techniques, one of which consists of studying the pipeline's response 
to signals introduced in the detector strain as software injections. This procedure helps 
understand the background noise present in the detector data and provides insight for 
fine-tuning the numerous parameters and thresholds in the different analyses. Since 
the strongest burst of gravitational radiation originates during the BBH merger, a 
good test to check the efficiency of existing DA techniques is to inject NR signals that 
contain the final part of the coalescence into simulated LIGO /Virgo data and proceed 
with the GW searches. 

An open collaboration under the acronym NINJA [17] has been created, that 
aims at fostering closer interaction between numerical relativists and data analysts. 
Its immediate purpose is to understand the sensitivity of current DA pipelines to 
binary black hole NR signals buried in simulated noise. To that effect, a set 
of waveforms was injected into LIGO- and Virgo-like noise, within the mass and 
distance ranges displayed in Figure [T] For the numerical waveforms that were used 
see [TBI US HOI US Ull HH HSl [HI H3 HSl I3 EOl EU, for descriptions of the 

numerical codes see[31[Ml[Ml[13[3i[Ml[li[371[3Hl[SSllini[SI]- Several data analysis 
groups took part in the NINJA project, employing their data-analysis codes to search 
for the NR signals. A comprehensive compilation of the goals and results arising from 
the NINJA collaboration can be found in a dedicated paper [41] . 

This paper presents the results of a search for BBH signals in the NINJA data, 
employing a new template bank based on a phenomenological model for non-spinning 
BBH systems by Ajith et al [T3J [JS] [TBJ . These waveforms fully describe the three 
stages of the coalescence of two non-spinning black holes for small mass ratios and, if 
implemented in current GW searches, should represent a step forward with respect to 
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PN templates that do not account for the merger and ringdown phases of the process. 
The remainder of the paper is organized as follows. In Section [2] we introduce the 
data-analysis pipeline developed by the Compact Binary Coalescence (CBC) group 
of the LSC to search for signals from compact binaries in the LIGO/ Virgo data. 
In Section [3] we describe the modifications performed in the pipeline that allow the 
use of the phenomenological waveforms as a new search template. In Section [4] we 
present the results obtained when applying this new template family to a search in 
the NINJA data. Finally in Section [5] we review our main conclusions and propose 
future improvements and directions for further research. 

2. The inspiral pipeline 

The LSC inspiral pipeline infrastructure developed by the CBC group has been 
employed to analyze the data released for the NINJA project. The pipeline, without 
major conceptual modifications, has been used in LSC searches for compact binaries 
from the third LIGO science run onward [lOi lllj . The same pipeline has been 
modified for the analysis of the NINJA data with the phenomenological template 
family described in Section [3] Since the signals present in the simulated noise are 
known to be NR simulations of BBH coalescences, the search method consists of a 
matched-filter technique [32] using an IMR waveform model based on the hybrid NR- 
PN waveforms as described below. The findings presented here concentrate solely on 
the simulated LIGO detectors HI, H2 and LI, although the NINJA data was generated 
for the simulated Virgo interferometer VI as well. 

The LSC inspiral pipeline performs a scries of hierarchical operations in order to 
search for real signals buried in the detector noise. The first stage is the construction 
of a suitable template bank covering the desired parameter space and for which the 
overlap between any point in the template space and the nearest template is at least 
as large as a certain minimal match [43l [44] , set as 0.99 in our search. After a filtering 
process through the desired template bank, the triggers that pass a certain signal-to- 
noise (SNR) ratio threshold are checked for coincidence in time and masses between 
two or three detectors. For our analyses the commonly-used threshold value of 5.5 
was employed. Whenever triggers are found with comparable coalescence time and 
parameters (in this case, component masses), they are stored as coincident [3S]. The 
detection statistic is the combined SNR of the single detector triggers. 

Once the initial matched filter has produced a list of triggers that pass the first 
coincidence step, a second stage follows where data is again filtered, but only through 
the templates that previously matched a trigger. Additionally the x 2 [IS] and r 2 [47] 
signal-based vetoes, designed to separate true inspiral signals from fluctuations in non- 
stationary noise, are applied. At this point an effective SNR p e s,i(pi,Xi) ^ s calculated 
combining the standard SNR with the x 2 value characterizing the mismatch between 
the spectral content of the template and the data. After further coincidence tests, 
the surviving triggers are listed as true gravitational wave candidates and constitute 
the output of the search. The significance of the triggers is based on the combined 
effective SNR, namely 
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Figure 1. Total mass and distance of the NINJA injected signals, with the colour 
code encoding the modulus of the dimensionless total spin \S\/m\ + S%/Tn%\ of 
the black holes. The total mass of the injected signals lies within the range 
36Mq < M < 46OM0 and they are located at a distance between 52 and 
480 Mpc. The vertical line delimits the mass parameter space with M < I6OM0 
that our phenomenological template bank covers. The full NINJA data set spans 
a duration of a little over 30 hours and contains total of 126 signals injected in 
simulated noise, 67 of which overlap with the parameter space of our template 
bank. 



and the triggers found by the pipeline allows for conclusions about the sensitivity of 
the analysis and the relative performance of the different template banks. 

3. The phenomenological template bank 

Most standard searches for gravitational waves from BBHs use the PN approximation 
of general relativity to construct banks of templates that account for the inspiral 
stage of the coalescence process, and the final ringdown can also be computed via 
perturbative techniques. However, the full calculation of the waveform in the merger 
stage requires numerical methods. These numerical simulations are in general rather 
expensive, and it is at present not feasible to model a coalescing binary over hundreds 
of orbits with sufficient accuracy. It is in fact also unnecessary to do so, because 
PN theory provides a valid description of the system when the black holes are 
sufficiently separated and the gravitational field is weak. Thus, a promising approach 
for constructing long waveform models covering the inspiral, merger and ringdown 
regimes is to stitch together the results of PN and NR calculations. 

One procedure for constructing such hybrid waveforms is presented in [TH fTSl US] , 
where PN and NR waveforms are matched in an appropriate regime (—750 < t/M < 
—550) prior to the merger (M is the total mass of the binary system in solar masses). 
Restricted 3.5PN waveforms at mass-quadrupole order are used for the inspiral phase, 
as given by equation (3.1) of [15]. For the numerical part, the model is based on 
long unequal-mass waveforms from simulations run by the Jena group using the BAM 
code [48] |32j [49]. These simulations span a range of mass ratios corresponding to 
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Figure 2. The figure shows two time-domain phenomenological waveforms from 
the template bank used in this search, corresponding to equal-mass binaries in 
the corners of our parameter space, namely (20 + 20)Mq and (80 + 80)M© BBH 
systems. The original and "whitened" 1501 waveforms are shown, with their 
amplitudes arbitrarily resized. The dotted and dashed vertical lines mark the 
points where the ISCO and light ring frequencies are reached. The LRD frequency 
is not shown, since it basically extends up to the full waveform. A matched-filter 
search that starts at 30 Hz and ends at the ISCO will not be able to pick the 
most massive binaries, since the inspiral phase of the coalescence falls below the 
LIGO interferometers' detection band. It is expected that the light ring and LRD 
frequencies, which extend up to the BH merger and ringdown respectively, will 
show improved performance at recovering high-mass signals. 



0.16 < r\ < 0.25, where 77 = [m\m,2) j M 2 is the so-called symmetric mass ratio of the 
binary system. The matching of PN and NR data is performed over an overlapping 
region, under the assumption that both approaches to the true BBH waveform are 
approximately correct at the late inspiral stage. Once the hybrid waveforms are 
constructed, they are fit to a phenomenological model determined entirely by the 
physical parameters of the binary system. This fit to an analytical expression is 
performed in the Fourier domain, assuming a functional dependence of the form 

u(f)=A oS (fy*° aU) , (2) 

with amplitude and phase given by equations (11) and (14) of [15]. Each waveform is 
parametrized by the physical parameters of the system, which in the non-spinning case 
are solely the masses mi and mi of the black holes. As a result of the matching and 
fitting procedures described above, a two-dimensional template family of waveforms 
that attempt to model the entire coalescence of non-spinning binary black hole systems 
has been obtained. 

The phenomenological template bank has been included in the LSC inspiral 
pipeline routinely used by the CBC group as a new waveform for filtering in the 
time domain. A search on the NINJA data has been performed, within a mass range 
of 20M Q < mi,m 2 < 80M Q for the component masses, with 40M o < M < 16OM 
for the total mass of the binary. A technical issue regarding construction of high mass 
waveforms prevented us from using a parameter space that completely overlaps with 
the mass range of the injected NINJA signals (see Fig. [1]), which might negatively 
affect the capability for recovering high-mass signals. Further searches for prospective 
NINJA projects will greatly benefit from inclusion of high-mass templates, and a 
fix for this issue is underway. The template bank is constructed using the standard 
second order post-Newtonian metric, and uses a hexagonal placement algorithm in 
mass space with a minimal match of 0.99 [44]. Note that we have not implemented 
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the new metric correspondent to the phenomenological waveforms, which might be 
significantly different than the standard 2PN metric; however, as we shall sec, the 
current template bank placement suffices for detection purposes and for moderately 
accurate parameter estimation. The number of signals that are recovered by the 
pipeline depends strongly on the choice for the upper frequency cutoff used in the 
matched filter integral, as we have observed in our investigations with the integration 
stopping at the ISCO (Innermost Stable Circular Orbit, r = 6M), light ring (the 
unstable circular orbit for photons orbiting a Schwarzschild black hole, r = 3M) and 
Lorentzian ringdown (lrdJI frequencies. 

In Figure [2] we show two waveforms from our phenomenological template bank, 
which correspond to equal-mass binaries in the corners of our parameter space, namely 
total mass M = 40M Q and 160M Q . Displayed are both the original time domain 
waveform and its "whitened" form [50], as the initial LIGO detector perceives it, 
and the relative amplitudes have been arbitrarily resized. The whitened waveform 
is computed as the inverse Fourier transform of the original signal multiplied by the 
function (/) in the frequency domain, where Sh(f) is the one-sided noise power 

spectral density of the simulated LIGO detectors. In each plot the vertical lines 
correspond to the ISCO and light ring frequencies. In our searches we have started 
filtering against the phenomenological templates at either 30 or 40 Hz and we have 
stopped the integration at the three frequencies discussed in the above paragraph. 
It is evident that whereas a cut at the ISCO frequency still retains a good portion 
of the inspiral signal for low-mass binaries, it is insufficient for higher masses. The 
light ring and LRD frequencies, on the other hand, extend roughly up to the BH 
merger and to the Lorentzian tail (from the decay of the quasi-normal modes of the 
ringdown), respectively, and are therefore expected to produce more efficient results 
for a matched- filter search for high-mass signals. 

4. Results 

4-1- Efficiency for detection 

The main results of our search for numerical relativity signals injected in simulated 
LIGO noise employing a phenomenological template bank are presented in Table [T] 
We show here a summary of the found triggers at different stages of the pipeline for 
several runs, with the starting frequency for the matched-filter integral being cither 
30 or 40 Hz and the integration stopping at three different frequencies -ISCO, light 
ring and LRD- displayed in ascending order. We have separated our results in two 
sections, according to performance in recovering the full set of 126 NINJA injections 
and the reduced set of 67 injections whose total mass falls below 160Af Q . This choice 
is motivated by the construction of the phenomenological bank discussed in Section [3] 
A time window of 120 ms has been used in order to cluster the triggers found by 
the pipeline in a single detector. Similarly the coincidence has been determined within 
a 80 ms injection window. Given these choices for the parameters used in clustering 
the triggers, we report recovery of 80/126 triggers in double or triple coincidence 
for the full injection set and 48/67 triggers for the reduced set with M < 160M Q . 
These are triggers that survive the second coincidence stage (including the signal-based 
vetoes) for our best run, which corresponds to the matched-filter integral starting at 

X This is 1.2 times the fundamental ringdown frequency of Berti, Cardoso and Will, Phys. Rev. D 
73 064030 (2006). 
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Template 


Phcnom 


Phcnom 


Phcnom 


Phcnom 


Frequency Cutoff 


ISCO 


LightRing 


LRD 


LRD 


Filter Start Frequency 


30 Hz 


30 Hz 


30 Hz 


40 Hz 


Complete set of 126 N1NJA Injections 


Found Single (HI, H2, LI) 


78, 54, 69 


94, 66, 90 


92, 61, 87 


93, 60, 86 


Found Coincidence 


59 


78 


81 


80 


Found Second Coincidence 


59 


80 


80 


79 


Reduced set of 67 NINJA Injections with M < 160A/ Q 


Found Single (HI, H2, LI) 


40, 17, 32 


55, 41, 50 


55, 41, 50 


56, 40, 50 


Found Coincidence 


30 


47 


48 


47 


Found Second Coincidence 


30 


47 


48 


47 



Table 1. Results of the search for NINJA signals using the phenomenological 
template bank. There were 126 injections performed into the analyzed data for 
HI, H2 and LI, 67 of which fell within the mass range of our phenomenological 
template bank (Af < 160Mq). We explicitly show that a much better efficiency in 
trigger recovery is achieved when the cutoff frequency is pushed beyond the ISCO 
frequency, up to the light ring and Lorentzian ringdown frequencies. Likewise we 
observe improved efficiency in finding the signals that lie within the mass range 
of our template bank. In both cases the signal-based vetoes have little influence 
in the rejection of triggers, confirming their efficiency in separating inspiral-like 
signals from other kind of glitches. 



30 Hz and ending at the LRD frequency. The number of recovered triggers in the full 
mass range is compatible with the results quoted by other participants in the NINJA 
project employing searches with higher-order corrections PN templates extended up 
to larger frequencies, such as the effective ringdown (ERD) and weighted ringdowryj 
(WRD) ending frequencies [41]. The efficiency of the search improves, however, when 
we restrict ourselves to signals with masses overlapping those of our template bank. 
It is worthwhile noting that among the triggers recovered by the pipeline we find 
not only non-spinning simulations but also signals with non-precessing spins, such as 
the CCATIE and BAMJ3BH (except the a/m = 0.25 case, see below) waveforms. This 
supports existing evidence for the fact that non-spinning templates should be able to 
detect non-precessing spinning signals with moderate individual spins. The capability 
of non-spinning templates for recovering signals with precession and large spins could 
however be compromised, as we discuss below. Due to the low statistics of the present 
analysis, these statements should be taken with the appropriate reservations. 

Figure [3] provides an overview of the found and missed injections corresponding to 
total mass below 160 Af©. The colour code encodes the modulus of the dimensionless 
total spin \S\/m\ + S^/m^l of the black holes, and gives an indication of the 
injections that significantly deviate from the non-spinning case modelled by the 
phenomenological waveforms. We observe how signals located at distances above 
350 Mpc are systematically lost, giving us an indication of the distance reach of the 
pipeline; nevertheless, several nearby injections are missed as well. In order to track 
down the missed injections in the mass region below 160M Q , a compilation of their 
relevant physical parameters and associated information is given in Table O For a 
description of the diverse NR simulations listed therein we refer the reader to [5T] . A 
similar analysis of the missed and found injections has been recently performed by the 
Birmingham group in |51j , applying Bayesian inference on the NINJA data using a 



§ An intermediate value between the ISCO and the ERD frequencies. 
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Figure 3. The figure shows found and missed injections in the mass region 
32Mq < M < 160Mq as a function of their total mass and distance for our best 
search, starting at 30 Hz and stopping at the LRD frequency. The circles represent 
triggers that were recorded as either double or triple coincidences after the second 
stage (including the signal-based vetoes), whereas the triangles represent missed 
injections. The colour code represents the modulus of the dimcnsionlcss total spin 
\Si/mf + S^/mjj of the black holes. 



Nested Sampling algorithm. The work of Aylott et al explores how different waveform 
families affect the confidence of detection of NR waveforms. Their Baycs factor B is a 
metric for assessing the level of confidence that a signal has been detected, and their 
defined thresholds for log 10 B allow for classification of the signals as found or missed. 
The IDs displayed in bold type in our Table [5] correspond to signals that are reported 
as missed by the Birmingham group in Table 2 of [5T]; 12 of our 19 missed injections 
are also lost by them, a correlation that seems worth following up. Future versions of 
the NINJA project will certainly benefit from combined searches and cross-checks of 
this kind between different DA methods. 

Signals with large eccentricity, such as the Princeton e > 0.5 PU_T52W run are 
invariably lost by our pipeline. Likewise, the phenomenological templates are not able 
to pick up signals with considerable spin, such as the equal-mass, spinning waveforms 
from RIT with individual spins a/m = 0.92 pointing along the z-axis and BAM_FAU 
with randomly aligned spins a/m = 0.75. Our template bank was developed to 
search for signals in which spin is unimportant and no precession is present, so these 
results are understandable. Further work targeted to incorporating spins within the 
phenomenological model is desirable and will be undertaken in the future. The missed 
GSFC signals with mass ratio 1 : 3 and 1 : 4 correspond to runs with few orbits before 
merger and moreover they are injected at total masses bordering on the edge of our 
template bank, which could explain them being lost. More bewildering is however 
the fact that the pipeline misses a couple of long equal-mass non-spinning simulations 
injected at rather close distances, such as CaltechCornell, which is also missed by the 
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ID 


NR 
Simulation 


Tot. mass 
(M ) 


Distance 
(Mpc) 


Eff. dist. 
H (Mpc) 


Eff. dist. 
L (Mpc) 


V 


Spin 

a/m 


Tot. Spin 

\s\ 


136 


RIT 


94.6 


444.5 


15831.9 


2941.3 


0.25 


0.92 


1.84 


141 


RIT 


75.6 


355.8 


1047.2 


746.6 


0.25 


0.92 


1.84 


142 


RIT 


129.7 


442.5 


2221.8 


1537.7 


0.25 


0.92 


1.84 


59 


CC 


69.7 


87.5 


469.2 


1573.0 


0.25 








41 


BAM_HHB_spp25 


112.0 


212.1 


1150.8 


802.2 


0.25 


0.25 


0.5 


47 


BAM_HHB_sppOO 


150.0 


131.0 


648.5 


908.0 


0.25 








27 


BAM_FAU 


77.8 


121.3 


764.1 


741.1 


0.25 


0.75 


1.43 


29 


BAM_FAU 


72.3 


185.1 


1325.7 


885.4 


0.25 


0.75 


1.43 


30 


BAM_FAU 


73.9 


351.1 


896.9 


771.4 


0.25 


0.75 


1.43 


114 


PU_T52W 


82.6 


133.1 


364.0 


320.2 


0.25 








116 


PU_T52W 


88.2 


164.2 


654.7 


533.7 


0.25 








118 


PU_T52W 


90.0 


86.0 


1055.4 


452.3 


0.25 








120 


PU_T52W 


108.6 


106.0 


202.0 


205.9 


0.25 








125 


PU_T52W 


96.8 


258.8 


463.8 


464.7 


0.25 








126 


PU_T52W 


105.8 


398.5 


1175.2 


1797.1 


0.25 








64 


GSFC_X4 


134.7 


149.6 


1320.4 


856.0 


0.16 








68 


GSFC_X3 


160.0 


220.8 


819.3 


1123.2 


0.1875 








76 


GSFC_X4 


158.0 


145.4 


722.2 


558.5 


0.16 








95 


Lean_cl38 


68.6 


246.5 


333.9 


407.7 


0.16 












Total number of 


missed injections: 19 









Table 2. Overview of the 19 missed injections with total mass below 160Mq 
for the best run reported. The ID column stores an index that identifies each 
of the injections of the NINJA set. The convention for the naming of the NR 
simulations can be found in 1411 . Among the missed signals we stress the presence 
of waveforms with eccentricity, large spins and precession and also those injected 
at distances further than 350 Mpc. Note that the boldfaced IDs correspond to 
signals also reported as missed in [51| . where a Bayesian inference search on the 
NINJA data using a Nested Sampling algorithm is presented. 



search reported in [51], and BAM_HHB_spp00. A look at the columns of Table[2]that list 
the effective distance in H and L (the distance to an equivalent source with optimal 
location and orientation) indicates that it might be the poor orientation of these 
injections that prevents the pipeline from finding them. Aside from these individual 
cases, which would need a careful follow-up that is below the scope of this paper, 
we can justify the rest of missed injections as those either placed at large distances, 
presenting large spin values and/or precession and containing few orbits before merger. 

Among the obvious improvements that a search with phenomenological templates 
could benefit from we can mention the following. Firstly, and once the technical 
issue with the generation of high-mass templates is resolved, the search would clearly 
improve with the use of a template bank that fully covers the parameter space of the 
signals searched for. Additionally, the inclusion of the fourth interferometer VI in our 
pipeline shall provide a larger number of recovered triggers, in the manner reported 
by the search using EOBNR templates that is described in Section 4.1.3 of [41]. Both 
improvements will be most likely incorporated to searches with the phenomenological 
template bank in future realizations of the NINJA project. 
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Figure 4. Figures showing the accuracy with which the chirp mass (left panel) 
and end time (right panel) of the reduced set of NINJA injections with total mass 
below 160Mq are recovered using the phenomenological template bank. In both 
plots, the colour scale is given by the total mass of the system. The chirp mass, a 
common quantity in data analysis which is defined as A4 = (mimap/^M) -1 ' 5 , 
is typically recovered within a 20% — 40% accuracy, depending on the chirp mass 
and total mass of the system. On the right panel we observe a few low-mass 
outliers that would need a more careful follow-up. Nevertheless, the results for 
parameter estimation with our IMR bank constitute a significant improvement 
over current LIGO/Virgo searches with standard PN templates. 



4-2. Efficiency for parameter estimation 

The number of found versus missed triggers is not the only relevant metric for assessing 
the performance of the standard GW searches. If astrophysically relevant statements 
are to be made from GW observations, the ability of accurately estimate the physical 
parameters of the measured signals is crucial. Figure 2] shows parameter estimation 
for the phenomenological search on the NINJA injections. The left panel displays 
the fractional difference for recovery of the chirp mass of the system, defined as 
M. = (mim2) 3 / 5 (M) -1 / 5 . The vertical colour bar encodes the total mass in solar 
masses. We report substantial improvement in parameter estimation with respect to 
LIGO/Virgo standard searches that make use of PN templates (see Figure 8 of the 
NINJA paper [H]), although it should not be forgotten that the PN searches reported 
in [5T] make use of banks with masses up to 90 Mq only. In our case the chirp mass is 
recovered within a 20% accuracy for values below 40M Q and ~ 40% for signals with 
larger chirp mass. The outlier that can be spotted at Ai ~ 30M Q corresponds to 
a CaltechCornell waveform injected at 132 Mpc with total mass 56.6M Q . For this 
particular injection the accuracy in parameter recovery is rather poor, and further 
work to understand this behaviour will be undertaken in the future. 

The panel in the right shows the accuracy in end time recovery of the found 
signals, with the colour code again displaying the total mass of the system. The sign 
convention for the At onc i corresponds to the injected minus the recovered parameters, 
so that a trigger that presents a positive value At cn d > indicates that the signal was 
really injected at a later time than the value recorded by our pipeline. Most of the 
signals are recovered at a time within a few hundredths of second from the injected 
end time value, with the outliers corresponding partially to signals with larger total 
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mass. A more careful study of the meaning of the end time when we are dealing with 
NR signals is needed in order to improve accuracy in estimation of this parameter. 

Even though the number of total recovered triggers for the phenomenological 
search on the NINJA data is similar to the results quoted by the standard PN 
searches, it is clear that the use of a full IMR template bank helps the estimation 
of the physical signal parameters. In view of these and other coincident results quoted 
in [41j [52], [51] , we conclude that searches that attempt to recover and estimate the 
physical parameters of BBH signals in the mass range 10 2 — 10 3 Mq would profit from 
using an IMR template bank that fully models the inspiral, merger and ringdown of 
the binary system. This is of crucial important for future LIGO/Virgo searches that 
aim at targeting coalescences of compact objects in the above-mentioned range, for 
which a full template bank adapted to arbitrarily high masses (and ideally also to 
non-zero spin values) needs to be developed. Attempts in this direction are already 
underway within the LIGO /Virgo collaboration. 

5. Conclusion 

We have presented the results from the first search using the standard inspiral pipeline 
modified to match-filter against the phenomenological template bank introduced by 
Ajith et al, an inspiral-merger-ringdown template that models the full coalescence of 
non-spinning black hole binaries of small mass ratios. This procedure has been directly 
applied to the search for numerically-simulated signals injected into simulated LIGO 
noise within the frame of the NINJA collaboration. We have tried several values for 
the cutoff frequency in the match-filter integral and found results that corroborate 
the need for pushing the integration to higher values than the ISCO if one wants 
to pick up signals that contain information about the merger and ringdown stages 
of the BBH coalescence process, which is very relevant for LIGO/Virgo high-mass 
searches. The total number of NINJA signals that we are able to recover using the 
phenomenological template bank is comparable to the results reported by other groups 
that participated in the NINJA project searching with PN templates extended to 
higher frequencies. Nevertheless, improved accuracy in parameter recovery is obtained 
with our phenomenological template bank. Moreover, it is expected that, had we 
constructed a bank covering the exact parameter space of the NINJA injections and 
included the fourth VI detector in our analysis, the number of recovered triggers would 
have significantly increased. This information will be used as starting point in future 
applications of the phenomenological waveforms to the search for inspiral signals. 

There is however a number of signals that are not found by our pipeline and that 
we have identified as injections falling outside the parameter space of our template 
bank, injections placed at distant locations above 350 Mpc, poorly oriented injections 
and injections corresponding to evolutions of black holes with significant eccentricity, 
spins and precession, plus a few outliers that undoubtedly need a more careful 
and dedicated study. More precise statements on found and missed signals would 
nevertheless require larger statistics in our data. 

There are several obvious improvements in the construction of the phenomenolog- 
ical waveform model that can be made. The phenomenological waveforms have so far 
been constructed only for the dominant mode of the comparable mass non-spinning 
case and it is not surprising that, for example, some waveforms with precession or 
eccentricity are not recovered. Searches for spinning objects are of the greatest as- 
trophysical importance, and phenomenological or other types of theoretical models 
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for these kind of systems are much needed. The phenomenological ansatz used to 
fit the merger portion (especially the power-law for the amplitude) needs to be ex- 
tended to higher mass ratios and spins. The error bars on the phenomenological 
parameters need to be better quantified, and the matching to post-Newtonian theory 
done as early in the inspiral phase as possible, and higher modes included. Notwith- 
standing, our results show a clear improvement in parameter estimation when the 
full inspiral-mcrgcr-rindown waveform is used as template. This result is confirmed 
by other reports within the frame of the NINJA collaboration. Future research on 
phenomenological models comprises development of adequate phenomenological wave- 
forms targeted at getting better detection and parameter estimation of signals from 
black-hole binaries in the current LIGO/Virgo data-analysis pipelines. Work on these 
improvements is in progress, and will be applied in future NINJA projects. 
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